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Abstract 

We provide a corrected Green's function for a polymer chain trapped in a two-dimensional 
anisotropic harmonic potential with a fixed boundary condition. This Green's function is a modified 
version of what Doi and Edwards first derived to describe the polymer chain confined in the tube- 
like domain of surrounding entangled polymers [J. Chem. Soc. Farad. Trans. II 74 (1978) 
1802]. In contradiction to the results found by lanniruberto and Marrucci (IM) when applying 
the Doi-Edwards Green function [J. Non-Newtonian Fluid Mech. 79 (1998) 225], we find that 
the stress-optical rule is violated for any tube potential either circular or elliptic if the corrected 
Green's function is used. The violation is due to the presence of the virtual springs to confine the 
chain in the tube rather than the anisotropy of the confinement potential. On the other hand, 
Doi and Edwards used their Green's function only for estimation of the monomer density along 
the primitive path where we find just a small correction. Since they did not use it for rheological 
calculations, the stress-optic rule appears to be safe for the Doi-Edwards model. 
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I. INTRODUCTION 



For polymer melts, it is well established experimentally that there is a proportionality 
relationship between the_stress tensor and the refractive index tensor as long as the polyrner 
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chains are Gaussian 



This linear relation is called the stress-optical rule (SOR) j^. 
For Gaussian chains, the traction force of a chain strand whose ends are constrained by 
entanglements is a linear function of the strand's end-to-end vector, so that the stress tensor 
of the polymer melt is proportional to the second-moment of the strand's orientation vector 
[sl . On the other hand, the refractive index tensor is proportional to the same second moment 
tensor, so that the SOR holds for Gaussian chains. If the polymer chain is stretched beyond 
the Gaussian regime, SOR is violated because of the breakdown of the linear relation between 
the .action and the e.d-tc.e„dle„gtK„ftKe*and a Q. 

The tube model of entangled polymers assumes that the lateral motion of the polymer 
chain in the melt is prohibited, as if the chains were confined in a tube-like domain j^. 
The origin of the stress of the material is usually assumed to be the traction force of the 
chain strand along the tube. But in 1998, lanniruberto and Marrucci (IM) considered the 
possibility that the pressure exerted by the chain strand on the confining tube wall niay 
contribute significantly to the stress tensor together with the traction along the tube [7|. 
They then investigated theoretically whether or not the SOR is obeyed in such a case. Based 
on the assumption that the equilibrium axial symmetry of the tube will be destroyed by the 
deformation of the melt, and the local geometry of the tube segment will become biaxial, 
IM considered the situation that the tube cross section is elliptical, or, in other words, the 
lateral pressure on the tube wall is anisotropic. Effects of the lateral pressure were taken into 
account by introducing "virtual springs" that connect all parts of the chain strand to the 
central axis of the tube. Then IM calculated the stress tensor and the polarizability tensor 
by applying a slightly modified Green's function that was first derived by Doi and Edwards 
(DE) to estimate the monomer density per length of the tube with circular cross-section js, 
App. A]. With the help of a "free confinement assumption" (more about this below) that 
each link of the chain does not suffer from the confinement potential, IM purportedly proved 
that the SOR is obeyed even when the contribution from the anisotropic lateral pressure (or 
confinement potential of the virtual springs) contributes to the stress. 

However, the Green functions given by DE and applied by IM are actually incorrect 



because these do not satisfy the boundary condition or continuity condition (see Sec. [IT|) . 
These Green functions can approximately describe an extreme case where the confinement 
potential is so large that the contribution of the chain traction along the tube cross-section 
is negligible compared to that of the confinement potential (see Sec. IIIip . But IM did not 
take this limitation into consideration in their attempt to prove the SOR jj]. 

In this paper, we derive the corrected Green function so that the boundary (continuity) 
condition is satisfied (in Sec. IIIip . and reconsider if the SOR is obeyed in IM's model of 
entangled polymer melts. We show in Sec. llVl that. if the polarizability tensor and the stress 
tensor are calculated based on the corrected Green function, and also if the calculation is 
appropriately conducted without employing the free confinement assumption, then the SOR 
is not obeyed due to the presence of the confinement potential — regardless of whether the 
confinement potential is isotropic or anisotropic (i.e., the cross-section of the tube is circular 
or ellipsoidal) and also irrespective of whether the confinement is strong or weak. The 
deviation of the stress tensor from the polarizability tensor is exactly the stress components 
originating from the virtual springs of the confinement potential. 

DE estimated the monomer density along the primitive path js], App. A] based on their 
Green function. Although it does not satisfy the boundary (continuity) condition, the DE 
prediction of the monomer density is valid as shown in Sec ilV C\ On the other hand, DE 
did not use the Green function for rheological calculations, so that the SOR appears to be 
safe in the DE model for rheology predictions. 

Use of the corrected Green function might also be of considerable importance in micro- 
scopic studies of the confining tube potential. There have been recent attempts to model 



the confining potential numerically 
dynamics 
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through molecular 



- Il3| . Some of these 



9|, find the potential numerica 
10|, or predict the potential from a more microscopic basis 
works have suggested that the potential might be anharmonic, in contradiction with what 
was assumed by Doi and Edwards and here. There is also interest in how the potential 

nn 

might be affected by flow p^, 115|. Therefore, it should be important to have the correct 
statistics arising from the simple harmonic assumption first to facilitate mapping from the 
atomistic level. 

In this paper, we use the same notation as in IM's paper Q except that the components 
of the vector R are denoted as {Rx, Ry, Rz) instead of {x,y,z), and the unit vectors along 
each axis are represented by {d^, Sy, 6^) instead of {x, y, z). 



II. DOI-EDWARDS' GREEN FUNCTION FOR CHAIN CONFINEMENT IN A 
HARMONIC POTENTIAL 



For the purpose of estimating the monomer density of an entangled polymer strand per 
tube length, DE described the confinement of the chain in the tube by the two-dimensional 
isotropic harmonic potential given by js], App. A] 



V=^w\Rl + Rl) (1) 

where ksT is the thermal energy, w represents the confinement strength, and -R is the 
position of a point (called the bead or monomer according to IM, though the chain is 
described as a continuum object) on the chain. Then the Green function G{R, Ro]n,nQ) 
satisfying 

d V ' 

dn 



G{R, Rq] n, no) = S{n — nQ)S{R — Rq) 



(2) 



6 dR^ k^T 

gives the statistical probability that the nth Kuhn step is near R, for a strand trapped in 
the tube by the potential ([T]), given that the noth bead is at Rq. The delta functions in the 
right-hand side take account of the boundary conditions G{R, Ro;n,no) = (for n ^ uq) 
and G{R, Rq] uq, no) = 6{R — Rq) P, p. 17-19]. It is important to note that the second delta 
function on the right-hand side of Eq. ([2]) makes the chain continuous 16|, Sec. 3]. b is the 
persistence length of the chain, and also represents the strength of the connector springs. 
We call the portion of the chain with the persistent length b the link, as did IM. DE used 
the Green function in Ref. sl, App. A] to show that the equilibrium monomer density per 
tube length is of the order of where a is the dimension of the tube cross- sect ion. We 
note that the conformation of the confined strand is governed by the Ornstein-Uhlenbeck 
process [l3| when n is interpreted as time because the confinement potential is harmonic. 
DE gave the following expression jsl, Eq. (A. 6)] as a solution to Eq. ([2]) B 

G<-'(fl,0;n.O) « ^exp + flj, - M _ !^^] . (3) 

But Eq. (|3]) does not satisfy the boundary condition G{R, 0; 0, 0) = S{R) for the x and y 
components of R, thereby violating the continuity of the chain. Moreover, in the absence of 
^ A factor 2 in the denominator of the first term of the exponential function is missing in Eq. (A. 6) of 



Ref. 8|, and a prefactor of the exponential function in Eq. (2) of Ref 



correction seems not to be reported in reference [14] of Ref. 
R, R' , n, n' , not just for R' = 0, n' = 0. 



7| is a typo. Also the misprint 



7[ . Note that the solution is used for arbitrary 
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the confinement potential {w = 0), Eq. (|3]) becomes uniform for x and y components while 
both ends are pinned (which is only possible because of the discontinuity of the chain). That 
is, while the ^-coordinates of the monomer are continuous, the other coordinates are not, 
so that monomers are free to fly apart along the direction perpendicular the tube central 
axis. (If w is large, they can stay relatively close, but the chain is still discontinuous.) These 
are results of the fact that the effects of traction along the x and y axes are not included 
in Eq. (|3]); only influences of the conflnement are taken into account. Thus Eq. (|3]) leads 
to the wrong stress and polarizability tensors that do not include the effects of traction 
perpendicular to the tube central axis. The derived two tensors are appropriate only when 
the conflnement is much stronger than the traction. This limitation should be kept in mind 
whenever one uses or applies Eq. ([3]). 

III. CORRECTED GREEN FUNCTION FOR ANISOTROPIC HARMONIC PO- 
TENTIAL 



Here we flnd the Green function that satisfles Eq. (jl]). Like IM, we consider the general, 
anisotropic chain conflnement whose potential is described by 



V='^(wlRl + w'Rl). 
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(4) 



Since the x, y, 2;-components are decoupled, the solution of Eq. ([2]) can be decomposed as 

G{R,Ro]n,no)= Y\_ Gp{Rp,Ro^p\n,nQ). (5) 

I3=x,y,z 

We show that the solutions for the perpendicular components /3 = x, y is (see the following 
subsection IIII A|) 

{n - no)hwp 



Gf}{Ri3, -Ro,/3; n, no) = J/j exp {^-^ coth 
where the prefactor is given by 



J. 



csch 



{n-no)bwi3 



2Tih 

and the flrst moment is 



exp ( J- coth 

26 



[Rh - Rfi)' 



[n - no)bwp 



(for (3 = x,y), (6) 



Rp = -Ro,/3sech 



{n - no)bw,3 



(8) 
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For the parallel component P = z, the solution is 

Gz{R,Ro,z]n,no) = 



exp 



3 — -Ro,; 



(9) 



One can see that Eq. ([6]) is a plausible solution of Eq. ([2]) for the following reasons. 
Firstly, in the weak confinement limit W/3 — (for P = x,y), Eq. ([6]) becomes the same form 
as Eq. ([9]), i.e., 



Gi3{R, Ro^f};n,no) 



exp 



3 {Rfs — Ri 



0,13) 



(10) 



2TT\n — nolb"^ \ 2\n — nolb"^ 

as expected for a Rouse chain [iS], §2.3]. Secondly, if n is close to uq for a fixed finite Wjs, 
Eq. ([6]) is also approximately written as Eq. ( ITOj) . Thus, in the limit of n — ?■ no, Eqs. ([6]) as 
well as dH]) become Si^Rp — -Ro,/?) (for f3 = x,y,z), thereby satisfying the boundary condition 
Gi3{R, -Ro,/?; '^o, ''^o) = ^(-R — Ro,i3) for all components. And thirdly, if the inequality 



nbwi3 ^ 1 (for (3 = x,y) 



(11) 



is satisfied, then Eq. (|5]) or 
G{R, 0;n,0) 



TT ^ 

/3=x,?; 



csch 



nbWj3 




3 


exp 



I3=x,y 



W0_ 

2b 



coth 



p2 _ 3-R2 



(12) 



becomes 



G(^^)(i?,0;n,0)oc^exp 



w. 



3Rl nb{w^ 



2y- 2f«»-ii-=^ (.6.,»1)(13) 

where, for simplicity, we put no = and .Ro = without loss of generality. Equation (fT3!) is 
the Green function that IM used (Eq. (2) of Ref. Q]) to derive the polarizability and stress 
tensors. For an isotropic potential {wx = Wy), Eq. (fT3l) reduces to Eq. ([3]). Thus the Green 
function of IM (or DE for the isotropic case) is appropriate only when condition (fTTj) is 
satisfied. 

The inequality (ITT!) has two interpretations. For a given strength bwjj of the confinement 
potential, Eq. ffTTl) indicates that n is so large as to reach the asymptotic state (or the 'steady 
state' of the Ornstein-Uhlenbeck process if n is interpreted as time) where the effect of the 
boundary condition is negligible. On the other hand, if n is given, Eq. ( fTTl) implies that 
the strength of the confinement potential is so strong as to overcome the traction for the x 
and y components. DE apparently neglected the traction along the x and y axes to obtain 
Eq. ©. 
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A. Derivation of Eq. ([6]) 



Gp is expected to be Gaussian because the potential is harmonic. Therefore we can 
assume that it has the expression 



Gfs = exp (^-f{n)Rl - g{n)Rp - h{n)^ . 



(14) 



Now we substitute this expression into the /3(= x, i/)-component of Eq. ([2]) but without the 
delta functions (which are enforced below) 

,2 



dn 6 



6 



(15) 



Equating like powers of i?^, we obtain three ordinary differential equations 



dh{n) 



-fin) 



-9[n) 



0. 



dn 3 '' ^ ' ' 6 
Equation (II Gal) is the Riccati equation for /(n) and has solution 



(16a) 
(16b) 
(16c) 



f(n) = -^coth 



{n - no)bwi3 



(17) 



where uq is the constant. It approaches f{n) 2(n-no)b^ unconstrained limit tW/j — )■ 

which is appropriate. We note that there is another solution of Eq. f ll6ap 

(n - no)bwi3 



fin) = — tanh 
■'^ ' 2b 



(18) 



However, this solution goes to at ^ 0, and is therefore unphysical since the delta- 
function boundary conditions are not satisfied. Therefore we employ Eq. (fT7|) as the solution 
of Eq. (]16ap . Then g{n) and h{n) are obtained from Eqs. ( ]16bl) and (]16cp as 

{n - no)bwp' 



g{n) = C^csch 



(19) 



and 



h{n) 



log ( sinh 



(n - no)bwp 



2wp 



coth 



[n - no)bwp 



(20) 



FIG. 1. A chain strand confined in the tube-hke region. Confinement potential due to the tube, 
described by the virtual springs in this figure, affects all portion of the segment in a continuous 
manner. Ends of the chain strand are pinned on the central z-axis of the tube segment. 



respectively, where C is a constant. Now that we have Gr, the first moment of Rr is 



calculated as 



JZRG,{R)dR b , 

•= n (TD\j-D = ~^ — ^^^^ 



{n - nQ)hwp 



Thus, by putting n = uq, the constant C is determined as 

^ — — -Ko,/?-^, 



(21) 



(22) 



and then the first moment is given by Eq. ([8]). Multiplying Eq. f lHj) by the constant -> , 
so that the boundary condition is satisfied, is written as Eq. ([6]). The ^-component 
(Eq. Qj) is obtained by replacing /3 with z in Eq. ([6]) and by taking the limit Wz — )■ 0. 



IV. APPLICATION OF THE CORRECT GREEN FUNCTION 

In this section, we calculate several quantities by applying the Green function given by 
Eq. ([6]), and examine the SOR. 

A. Number of monomers in the tube 



We consider the simple case that the start and end of the chain strand are pinned at 
the central axis of the tube separated by the distance /, as both DE and IM assumed (see 
Fig. [1]). By putting Rq = and R = 18 ^ in Eq. (|6]) (or substituting R^ = Ry = 0, Rz = I 
in Eq. ( fT2|) ). we have 



G{R= /5^,0,n,0) = exp 



I3=x,y 



2'Kh 



csch 



nhwr- 



3/2 1, / 3 
h - log 



.(23) 



The most probable value of n is obtained by maximizing Eq. (123|) with respect to n for a 



fixed I 



a 



If we take the derivative of the argument of the exponential function with respect 



to n and set the result to 0, we obtain 



n = (24) 



b\jbT.li=x,y coth 



nbwfj 



In the condition that the inequality f|TT]) is satisfied, coth.[nbWj3/3] is approximately 1, and 
3/n in the denominator can be neglected. Then the most probable value of n that DE j^, 
Eq. (A. 9)] and IM [3, Eq. (3)] derived for the isotropic and anisotropic potentials, respec- 
tively, is obtained, i.e., 

3/ 

n ~ — _ (for nbwj3 ^ 1) (25) 

bV2bw 

where we put w := {w^ + Wy)/2 as the mean strength of the confinement potential. On the 
other hand, in the weak confinement limit, Eq. ( |24|) becomes 

n ^ ^ {w^^ 0). (26) 

An expression that interpolates Eqs. (!25|) and (126|) gives an approximate solution to Eq. 

as 

n ^ , (27) 

by/2bw + 967/2 

See Fig. |2]for reference. By comparing Eqs. (125|) and (126|) . the crossover of w between the 
weak and strong confinement can be estimated as U7c — If the condition Wj3 <^ Wc is 
satisfied, then the inequality opposite of Eq. (fTTl) is fulfilled because nfety^ ^ nb'^/P ~ 1. 



B. Tube dimension 

The conditional probability density of the kth bead's position Rk when the strand's head 
and tail are pinned at Rq = and R = ISz, respectively, is given by the product of two 
Green functions as 

Pk{Rk\Ro, R) oc G{Rk, Ro] k, 0)G{R, Rk] n, k). (28) 
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FIG. 2. The equilibrium number of monomers of a chain strand in the tube plotted against 
the strength of the confinement potential. For simplicity an isotropic potential Wx = Wy{='- w) is 
considered. Solid line is the exact result (numerical solution of Eq. (j24p ). dotted line is the IM (or 
DE) result (Eq. (p5|) ). and dashed line is the approximate result (Eq. (p7|) ). 




FIG. 3. Squared radius of gyration of the tube cross section along the axis perpendicular to the 
tube plotted against the strength of the confinement potential. For simplicity an isotropic potential 
is considered (a^ = ay =: a). Solid line is the exact result (Eq. (j3ip ). dotted line is the IM (or DE) 
result (Eq. (f32l) ). and dashed line is the approximate result (Eq. (f33|) ). 



The second moment of the k th bead's position along the x-axis is calculated as 



{Rk,x) 



I ^k,x Pk{Rk\Ro, R)dRk 

J Pk{Rk\Ro, R)dRk 
b 



Wx (coth [^] + coth 



{n—k)bwx 



(29) 
(30) 
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The radius of gyration of the elhptical cross-section of the tube along the x axis can be 
estimated by taking the average of the last equation over all k, i.e., 

where L{x) := coth(x) — 1/x is the Langevin function. A similar relation holds for the 
y-component. If inequality ( fTTl) is satisfied, then Eq. ( I3T1) is written approximately as 



al ~ (for nbw^ > 1). (32) 



This is the result that DE (for the isotropic case) [8|, Eq. (A. 7)] and IM (for the anisotropic 
case) 0, Eq. (5)] derived. Equation (|32l) corresponds to the second moment for the Ornstein- 
Uhlenbeck process in the steady state (when n is interpreted as time). In the absence of 
the confinement potential {w^ = 0), Eq. (13 ip reduces to = /^/18 as it should. By putting 
Eq. (p7|) into Eq. (!3T]) . we obtain the approximate expression for the tube dimension as 

^ ). (33) 

See Fig. [3] where these approximate results are compared with the exact one. 



C. Monomer density 

DE showed that the monomer density per tube length, n/l, is of the order of by 
using the Green function ([3]) which is valid only when the confinement is strong . That is, 
by eliminating w from two relations n ~ 3l/{bV2bw) (Eq.(l25])) and a ~ ^b/{2w) (Eq.(|32l)) 
for the isotropic case, the monomer density can be estimated as n/l ~ 3a/b^. Interestingly, 
this prediction is valid even when the confinement is weak. By using two relations n ^ 
(see Fig. [2]) and a ~ //(3-\/2) (see Fig. [3]) that hold in the weak confinement case, we obtain 
a relation n/l ~ 3\/2a/b'^. The difference in these predictions is just a prefactor of order 
unity ("\/2), and we can therefore conclude that the DE prediction for the monomer density 
is appropriate for any strength of the confinement potential (see Fig. H]). 



D. Expression of the polarizability tensor 



Here we derive the polarizability tensor A from the tube segment. In the calculation of 
the polarizability tensor, IM assumed that each link of the strand is not under the influence 
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DE's result 






approximation 





















lo--* 10-^ 1 10^ 

potential strength, wfi/b 
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FIG. 4. Monomer density (scaled by a/b'^) as a function of the strength of the confinement 
potential. 

of the confinement potential (what we call the "free confinement assumption"). On the 
other hand they did not use this free confinement assumption in the calculation of the stress 
tensor. Here we avoid the free confinement assumption for the calculation of both tensors 
because it seems inconsistent with Eq. ([2]) and is anyway unnecessary. 

The origin of the polarizability and the stress of polymer melts is the orientation distri- 
bution of the end-to-end vector of units of the chain. Like IM, we assume that this unit is 
each link (i.e., a segment of length b along the continuum chain). Validity of this assumption 
is discussed at the end of this section. Then the local polarizability tensor from the k th 
individual link is given by 

Ak = c{{Rk+i - Rk){Rk+i - Rk)) (34) 

where c is a constant. The brackets (■ ■ ■ ) indicate taking the average for the conditional 
probability density that the k th link's head (tail) is located at Rk (Rk+i) when the ends 
of the chain strand are fixed at Rq = and -R = Idz as before (see Fig. [1]) 

Pk{Rk, Rk+i\Ro, R) oc G{Rk, Ro; k, 0)G{Rk+u Rk] k + 1, k)G{R, Rk+i, n,k + 1). (35) 



By using this conditional probability density, the xx-component of Ak is calculated as 

/ dRk J dRk+i {Rk+i,x — Rk,xY Pk{Rk, Rk+i\Ro, R) 



k+l,x 



Rk,x) ) 



J dRk J dRk+i Pk{Rk, Rk+i\Ro-, R) 



(36) 



— csch 


nbwx 


sinh 


bwx 


sinh 


kbwx 




3 




3 




3 



sinh 



[k + 1 — n)bwx 



X coth 



{k 



n)bwx 





kbwx 




bWx 


— coth 




— 2tanh 




3 




6 



(37) 
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The i/y-component satisfies a similar equation. 

Now we calculate the principle components a^x, ttj/y, «zz of the polarizability tensor A 
from the entire chain-strand confined in the tube segmentn. Since n ^ 1, this is given by 
integrating A^. over k for all links of the chain strand. The x-component can be obtained 
from Eq. f l37p as 

: c j {{Rk+i,, - Rk,xf)dk = c— (38) 



a.7 



where 



6 



bWr, 



sinh 



6 



csch 



nbwr 



sinh 



(2n — l)bw^ 
6 



nbwx 



cosh 



bw^ 



sinh 



bw^ 
6 



The summation of Aj. over k yields a similar result but without the second term in the 
parentheses of [fj. A comparison of this expression of with IM's result is discussed in 
App. |Al The y-component has a similar expression. On the other hand, the ^-component is 

the polarizability from the Gaussian strand free from the confinement potential, i.e., 

nb"^ 



a. 



(42) 



where 



, 1 3/2 
1 \ ~ 

n 71%"^ 



1 + 



3/2 



(43) 



As a result, the polarizability tensor of a single tube segment aligned along the z-axis is 



(44) 



P=x,y,z 

IM took an average of the polarizability tensor A for a single tube segment over all tube- 
segment orientation and length, and also multiplied it by the number density of the tube 
segment in the melt to obtain the polarizability of the material. They did the same for the 
stress tensor. In this paper, however, we don't perform these procedures for both quantities 
because, for the purpose of confirming the SOR, it is enough to compare the polarizability 
^ In this paper we arc following the IM assumption that the stress and polarizability tensors originating 
from the chain strand in the tube segment have the principle directions, one of which is along the central 
axis of the tube segment (z-axis) and the other two {x and y axes) are along the principle directions of 
the ellipse cross section. 



n-l 

O-xx = c'^{{Rk+l,x 

k=0 



Rk,x) ) =c — 1 - 1/n + csch 



nbwx 


sinh 


{n - l)bwx 


^ tanh 








3 


3 




6 



2nb . ^ 
c sinh 



bwx 



csch 



nbwr 



sinh 



(2n - l)bwx 
6 



(40) 
(41) 



where 1/n in the parenthesis was neglected to obtain the final result. 

13 



and stress tensors only for a single tube segment because the origin is the same. Also IM 
replaced n in the expression of the polarizability tensor with its most probable value, which 
is, in our case, given as a solution of Eq. 

Finally we comment briefly on the proper choice of the smallest unit responsible for the 
polarizability and the stress of the material. IM assumed that h is the unit of polarizability, 
but here we consider the infinitesimal unit in accordance with the current, continuum picture 
of the chain. If we had discretized the chain by A as the unit of the polarizability, then the 
local polarizability from this portion would be given by A^'^^ := c{{R(^k+i)A—RkA){R{k+i)A — 
RkA))- Integrating A^^^ over all these portions < A; < n/A, the total polarizability tensor 



is obtained as A^^^ = c 



n 



(A) 



hwp/S. 



sinh 



3 E 
hwR^ 



(A) 
x,y,z ^^f) 



U},^' dpSfi where 



6 



X ^csch 


nhwj3 


sinh 


'{2n - A)hwp 




3 




6 



nhw^ 



cosh 



bwxA 


sinh 


bw^A 


3 




6 



1-^ 

n 



3/2 



:A 



31' 



-A. 



(45) 
(46) 



The previous result (Eq. f l44p ) is recovered by putting A = 1. On the other hand, in the 
limit of vanishing discretization A — for a fixed finite b (or the continuum limit), the 
polarizability matrix becomes isotropic 

nb"^ 



lim A(^) 



c S, 

3 ' 



(47) 



where 6 is the unit tensor. The isotropy comes from the elimination of the wavelength cutoff 
which makes the contour length of the strand infinite, as originally modeled. Since we are 
not interested in this trivial result, we discretized the otherwise-continuous chain by putting 
A = 1 as did IM. 



E. Expression of the stress tensor 

Here we calculate the stress tensor, and compare the result with the polarizability tensor 
derived above. 

IM discussed that the contribution to the stress is two-fold; one is the traction of the 
chain strand and the other is the confinement potential. But they considered the traction 
only along the tube segment (i.e., 2;-axis), which is consistent with their Green function. 
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Here we take the traction of all directions into account to be consistent with our Green 
function. 

According to the conventional theory for Gaussian chains 5|, |6|, the stress tensor arising 
from the chain traction is given as the integral of the second-moment tensor of the link's 
end-end vector over k. The contribution from a single tube segment aligned along the 2;-axis 
is given by 

T, = nk^T ^pSp^P (48) 

where 11/3 is given by Eqs. f l39|) and f H3|) . Therefore Tc is proportional to A given by Eq. (jH]), 
indicating that the SOR holds if the stress does not include the contribution from the virtual 
springs representing the confinement potential (or 'intrachain pressure'). 

On the other hand, the force arising from the virtual spring along the /3 (= x, ?/)-axis is 
■^1^ ~ ~ regarded this force as the origin of the pressure on the tube 

wall. Thus the contribution to the stress tensor from the virtual springs is 

V / {Fk,^Rk,p)dk S^dp (49) 
E I (^m)^^ ^^^^ (50) 



=x,y 

k^Tw^ 

P=x,y 

X ^ nl 

~6 "V 3 



= -I-bT Y. (51) 

I3=x,y 

where we used Eq. 0311) . Consequently, by adding these two parts, the total stress tensor for 
the tube segment is obtained as 



Te + T, 

" ' '"^-Lri^^lv. + nA4| (52) 




6 

In the process of deriving Eq. ( 152|) . one sees that the stress tensor T is not proportional 
to the polarizability tensor A because of the presence of T„ in the stress. The virtual spring 
contributes directly to the stress, but not to the polarizability; it affects the polarizability 
indirectly only through the Green function. Thus the stress-optical rule is not generally 
obeyed in the present theoretical model of entangled polymer melts by IM. It should be 
noted that the violation of SOR is not caused by the anisotropy of the potential but by 
the presence of the potential itself. That is, the SOR is violated even if the confinement 
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potential is isotropic [wx = Wy). Also, one might expect from the IM prediction that the 
SOR holds if the confinement is very strong because IM's prediction is based on the Green 
function that is appropriate in the strong confinement case. However the SOR is violated 
even when the confinement is strong due to the presence of the confinement potential as 
shown in App. [Bl 

We make a final note about the interpretation of the wall pressure. The sign of was 
chosen to be the opposite of so that is an outward force (i.e., 'pressure') toward the 
tube wall, versus Tc representing the traction. However it is unclear to us how to describe 
such a pressure by the confinement potential Eq. (j4]) which is, by definition, supposed to 
describe the attractive force toward the central axis of the tube. We do not go into details 
of this issue, but rather just conclude this section by remarking that the SOR is not obeyed 
irrespective of the sign of . 



V. CONCLUSIONS 



We derived the Green function for an entangled polymer chain trapped in a tube having 
ellipsoidal cross-section described by an anisotropic harmonic potential. Unlike the Green 
function derived by Doi-Edwards and that modified by lanniruberto-Marrucci (IM), ours 
satisfies the boundary condition along the axes perpendicular to the tube central axis. The 
stress tensor and polarizability tensor derived from our Green function without the free 
confinement assumption do not satisfy the stress-optical rule in the model proposed by IM. 
The stress-optical rule is violated because the virtual springs of the confinement potential 
contribute only to the mechanical stress tensor, not to the optical polarizability. Thus the 
presence of the virtual spring itself, rather than the anisotropy of the spring potential, is 
the source of the breakdown of the stress-optical rule. 



We are grateful to the Army Research Office (grants W911NF-08- 2-0058, W911NF-09-1- 
0378 and W911NF-11-2-0018) for financial support. 



16 



Appendix A: Comparison of IM's polarizability and our expression 



With the help of the free confinement assumption, IM derived the polarizabihty from 
their Green function which is appropriate only for nhw^ ^ 1 as [3, Eq. (17)] 

On the other hand, our exact expression (Eq. ( !38l) or Eq. (l4Tl) ) can be written approximately 
for the same condition as 



(^xx — C 



!^ (1 _ 6-^-^/3) (for l/n < 6w J . (A2) 



(It should be noted that Eq. ( 1A2I) divided by n corresponds to the mean-square displacement 
of the particle of the Ornstein-Uhlenbeck process in 'steady state' with the 'lag time' cor- 



responding to a single link 



16 



p. 77].) It is natural that Eq.(jAT]) and Eq.(lA2]) are different 
because these were derived under different assumptions. But surprisingly, both equations 
are equal up to first order in bw^ 

aif ) = a,, ^ (l - ^) (for l/n « bw, « 1). (A3) 

This coincidence is due to the cancellation of errors in axx^^ from (i) their Green function 
that gives rise to errors for small bw^ and (ii) the free confinement assumption. The higher- 
order terms of a^^ and a^i^^ do not agree because of the free confinement assumption that 
makes the Wa;-dependence weaker. See Fig. |5] where Eqs. (IA1|) . (IA2p and (IA3P are compared. 
A similar discussion holds true for the y-component. 

Indeed Eq. flASp happens to hold in the condition l/n bw^ <^ 1, but IM compared 
this result with the stress derived under the different conditions {l/n ^ bw^ and bw^ ^ 1) 
in order to confirm the SOR. If these quantities are compared at the same conditions, we 
see that the SOR is violated as shown in AppjBl 



Appendix B: Investigation of IM's discussion on the stress-optical rule 

In this appendix, we consider in detail how we arrive at the conclusion that conflicts with 
IM. For this purpose, we focus on the 'steady state' (1 ^ nbwp) as IM implicitly did. For 
clarity, we consider three cases where the confinement strength is (i) just lower bounded 
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FIG. 5. Polarizability plotted against the strength of the confinement potential. Solid line is 
Eq. ()A2p . dotted line is the result IM obtained, and the dashed line is the approximate curve for 
both lines (Eq. dMl)). 

1/n <^ hwp < oo, (ii) relatively weak \/n <^ hwp ^ 1, and (iii) very strong 1/n ^ hwp, 
1 < hwp. 

(i) If 1/n <^ hwp < oo, Eq. fl52l) can be approximated as 



T ~ nk^iT ^ 

\l3=x,y 



hWr- 



(1 



~bwfj/3 



) 



bWf3 

6 



SsSf 



no 



(Bl) 



and also Eq. (jUj) can be written as 



A c- 



\l3=x,y P ^ 



3/2 



(B2) 



where we used Eq. flA2p . SOR is violated due to the contribution to stress from the virtual 
springs (underlined term in Eq. ( IBll) ). 

(ii) \il/n <^ hwp <^ 1, Eqs. fIBip and (lB2p can be expanded in a Taylor series in bwi3, 
and consequently these can be decomposed into the isotopic and the anisotropic parts as 

S/2 \ 

(B3) 



bw,. . J. bw,. 

^xOx 



and 



A-^c—6 + c—(-^5 6 -^-^66 



3/2 



(B4) 



3 ' ■ ' 3 V 6 '■''-^ 6 '''' ' 
respectively. As discussed in App. \^ Eq. (lB4p for a single tube segment happens to corre- 
spond to the polarizability tensor that IM derived 0, Eq. (21)]. The anisotropic components 
of Eqs. (IBSp and (IB4p are not proportional to each other, thereby violating the SOR. 
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(iii) If 1/n <^ hwp and 1 ^ hwp, then the confinement is much stronger than the traction. 
Therefore the contribution from the chain traction — exp[— 6w^/3]) can be dropped 

from Eq. (IBll) . and the first term n in the parenthesis for the z-component can be neglected 
compared with the x and y components in Eq. fIBip . Thus we have 



This corresponds to the resuh that IM derived [Tl, Eq. (9)]. IM compared Eq. (1B5I) to the 
anisotropic part of Eq. (]B4p to conclude that the SOR is obeyed even in the presence of the 
anisotropic confinement potential. However, since the range of where each equation is 
approximately correct is not the same, one cannot verify SOR from such a comparison. The 
SOR is not satisfied within the present condition because the polarizability tensor has only 
the z-component as 



Thus we conclude that the SOR is not obeyed in this model except the case where the 
potential is so weak as to be negligible compared with the chain traction, and consequently 
both tensors are nearly isotropic: T oc A oc (5. But we are not interested in this extreme 
case because the chain strand is approximated by a free Gaussian chain without restriction 
as shown in Fig. [2] and Fig. [3], and therefore there is no surprise that the SOR holds. 
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